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Sur les Méthodes d'Approximations Successives dans la 
Théorie des Équations Différentielles.* 



Pae Emile Picard. 



J'ai consacré plusieurs Mémoires à l'application de méthodes d'approxima- 
tions successives pour démontrer l'existence et faire la recherche des intégrales 
de certaines équations différentielles, quand des conditions aux limites sont 
données qui définissent ces intégrales. Ces méthodes s'appliquent aux équations 
différentielles ordinaires comme aux équations aux dérivées partielles, mais pour 
ces dernières les conditions d'application sont bien différentes suivant que les 
équations considérées appartiennent au type elliptique ou au type hyperbolique. 
Les premières se rencontrent surtout en Physique mathématique et dans la 
théorie des fonctions ; je ne m'en occuperai pas dans cette Note.f Relativement 
aux équations du type hyperbolique, l'utilité de ces méthodes est d'une double 
nature. Elles permettent d'abord de faire la recherche des intégrales en suppo- 
sant les équations différentielles définies seulement pour les valeurs réelles des 



*Nou8 reproduisons ici une note insérée par M. Picard à la fin du tome IV de la Théorie des Surfaces 
de M. Darboux (Note de la Rédaction). 

t Relativement aux théorèmes généraux relatifs à ce cas, nous énoncerons seulement la proposi- 
tion suivante (Journal de l'École Polytechnique, 1890). Soit l'équation linéaire 

dx^ dxBy Sy^ dx ay ' 

où les coeflaoients sont des fonctions analytiques des deux variables réelles xeiy: toute intégrale de cette 

équation bien déterminée et continue ainsi que ses dérivées partielles des deux premiers ordres dans une 

région du plan pour laquelle 

B2— AO<0 

est une fonction analytique de x et y. Il est clair qu'il peut en être autrement dans une région où 
B^ — AC serait positif. On trouvera un théorème plus général (Comptes Rendus, Juillet 1895). 
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variables, et en faisant ainsi le minimum d'hypothèses sur ces équations ; c'est là 
un point d'un certain intérêt philosophique. 

Une conséquence pratique en découle ; on obtient, en général, pour les inté- 
grales un champ de détermination plus étendu qu'avec les méthodes fondées sur 
l'emploi de fonctions majorantes quand ces méthodes sont applicables. 

1. Rappelons d'abord, sans y insister, les résultats i-elatifs à une équation 
ordinaire du premier ordre 

Si/(a;, y) est une fonction réelle et continue des deux variables réelles 
ce et ^ , quand celles-ci varient respectivement dans les intervalles 

(»o — a, ajo + a), (.?/o — ^ . «/o + &) , 

et si, de plus, il existe une constante positive k telle que 

\f{x, y + Ay) —f{x,y)\<h\ Ay \ , 

x,yety-\- Ay étant les intervalles indiqués, et qu'enfin M désigne le maximum 
de la valeur absolue de /{x, y) dans ces mêmes intervalles, les approximations 
successives donnent l'intégrale de l'équation prenant pour x = Xo la valeur î/o, 
sous forme de série convergente dans l'intervalle {xq — p, Xo + p), en désignant 
par p la plus petite des deux quantités * 

a et -^ . (1) 

M. B. Lindelôf, qui a très heureusement approfondi cette question (Journal de 
Math., 1894), a même indiqué un autre champ de convergence qui peut quel- 
quefois être plus étendu que le précédent. Désignons par Mq la plus grande 
valeur absolue de f{x, y^, quand x varie de Xq — a à a^o + a ; un champ de con- 

*Nous avons supposé la fonction /(«, 2/) définie seulement pour les valeurs réelles de xety. Dans 
le cas où /(a;, y) est une fonction analytique de a; et 2/, holomorphe dans les cercles de rayons a et h 
tracés respectivement autour des points x^ et yo , et en désignant par M le module maximum de / dans 
ces cercles, les approximations successives permettent d'obtenir l'intégrale prenant pour x=.Xa la 
valeur y„ sous forme de série convergente dans un cercle de rayon f autour de «o (en désignant par p 
la même quantité que ci-dessus). La méthode des fonctions majorantes donne un champ de conver- 
gence moins étendu. 
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vergence assurée est l'intervalle {xq — p', x^ + p') , en désignant par p' la plus 
petite des deux quantités 



a et ^ 
k 



F'°<'+f)' <^) 



et p' peut dans certains cas dépasser p. 

Il n'est pas sans intérêt de rappeler que la première méthode de Cauchy, 
telle que nous la connaissons par les lepons qu'a rédigées M. Moigno, et qui a 
été depuis reprise par M, Lipschitz, méthode dont le principe est de considérer 
l'équation différentielle comme une équation aux différences, définissait précisé- 
ment l'intégrale dans l'intervalle correspondant à (1). 

2. Considérons maintenant une équation aux dérivées partielles de la forme 
d^z Tj, / dz dz 



F 



/ àz dz \ 



dxdy \ ' dx ' dy 

Les approximations successives permettent, entre autres problèmes, de 
former l'intégrale d'une telle équation se réduisant, pour œ = Kq, à une fonction 
donnée de y, et pour y=zy^ h. une fonction donnée de x. Je prendrai d'abord 
le cas de l'équation linéaire 

d^z dz . T dz . 

= a 3 1- b -g— + cz , 



dx dy dx dy 

où a, h, G sont des fonctions des deux variables réelles ce et y. Nous les sup- 
poserons continues à l'intérieur et sur le périmètre d'un rectangle R de côtés a 
et /? parallèles aux axes et dont (a-o, 2/o) sera le sommet de moindres abscisse et 
ordonnée. On veut trouver l'intégrale de cette équation se réduisant pour 
y=L y^ à 4> (as) et pour a; = «o à 4' («/) • La fonction ^ (œ) est continue de a^o à 
cco + a, et i^ (y) est continue de 2/0^2/0 + /?; on a, bien entendu, ^ {x^ = 4 (yo) 
et les deux fonctions 4» et i^ ont des dérivées premières continues. 
Envisageons, en premier lieu, l'équation 

oii/(a;, y) est une fonction donnée. La fonction 

12 



90 Picard : Sur les Méthodes d'Approximations Successives 

est l'intégrale de cette équation s'annulant, pour x =: ccq, quel que soit y, et pour 
2/ = 2/0 quel que soit x . Ceci posé, nous formons les équations successives 



dxdy 



= 0, 



2 =a_-i- -f 5» -}-c3i, 



3a; 3?/ 3a3 dy 



3a; 3^^ dx dy 

On intégrera la première équation en cherchant son intégrale % se réduisant 
à ^ (ce) pour î/ = 2^0 et à 4 (v/) pour as = ccq, intégrale qui est visiblement 

2i = 4» (*) + 4- (2/) — 4> (a^o) • 

Pour toutes les autres fonctions 2»(« >» 1), elles sont supposées se réduire à 
zéro pour as =: aso quel que soit y , et pour y = yo quel que soit x . 
Nous allons montrer dans un moment que les séries 

dzi , 3z2 , , 32„ ■ 

3^^ 3^ ^•••• + '3^ + ---" 

321 , 32g , , 3z„ 

sont uniformément convergentes dans le rectangle R. Ce point admis, on voit sans 
peine que la fonction 

Z = Zi + sg +....+ z„ + ... , 

est Vintégrale cherchée. On tire, en effet, des équations précédentes 

+ 5 gi ^i + • - • + -^^ + c (21 + . . . . + g« - 1)] ^x tfy , 

d'où l'on conclut à la limite, en s'appuyant sur la convergence uniforme des 
séries écrites plus haut, 

Z=::^{x) + ^{y)-^{x,) +£J\a^^ + hf^ -^ cZ]dxdy 



et enfin d^Z 3Z , , 3Z , „ 

oa; dz/ ox oy 
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Abordons donc la question de convergence. Je désigne par M le maximum de 

dans R, et par à; le module maximum de a,b, c dans ce même rectangle, et je 
considère le système 



dxdy 



M, 



ox oy ax oy 

9^«« _ t. 9^»-i 4. T, ^^n-l , 7.,, 

ax oy ox oy 

tous les u s' annulant pour x=^Xq quel que soit y, et pour y^^y^ quel que soit x. 
Si nous prouvons la convergence uniforme de la série 

Ml + «2 4- • • • • + W» + . (3) 

la convergence de la série des s en résultera immédiatement, car |2„|<C |'M„_i|. 
Or, soit 

on aura 9^Ui 



dxdy 



:M, 



1 = ^1 +^1 +Ui. 



dx dy dx dy 

dxdy dx dy 

Si la série des u est convergente, la série 

Ui + ^U3+.... +^''U„ + i+ .... (4) 

représentera l'intégrale de l'équation 

^'u -/fc|^ + >fc|^ + ;fcu + M, (5) 



dx dy dx dy 

s'annulant pour x = Xq, quel que soit y, et pour y=-y^, quel que soit x. Or si 
nous montrons que, pour l'équation précédente, l'intégrale satisfaisant à ces con- 
ditions initiales, est une fonction holomorphe de h pour toute valeur de h, la 
convergence de la série (3) sera établie, car cette intégrale devra nécessairement 
avoir la forme (4). 



92 Picard : Sur les Méthodes d'Approximations Successives 

Or l'équation (5) est facile à discuter. Prenant a-o = yo== 0> "^ows, poserons 

L'équation (5) devient 

J!|- = (P + A;)V + Me-*<»' + ^), (6) 

àxoy 

L'application des approximations successives à cette dernière équation est 
immédiate. On a à considérer les équations 

^ * 3^ Me"*" '" + ^' 
dxdy 

tous les V s'annulant pour a; = , ainsi que pour y = 0. En désignant par N la 
valeur absolue maxima de Vo dans R, on aura 






(1.2 nf 

d'oii l'on déduit de suite la convergence de la série 

V„ + Vi +.... + V„ + .... , 

qui représente l'intégrale cherchée de l'équation (6). La méthode des approxima- 
tions successives donne donc pour l'équation (6) une série convergente, quand 
(x, y) est dans R. Chacun des termes de cette série est une fonction holomorphe 
de ^, et la série converge uniformément, quel que soit k, dans un domaine fini 
quelconque du plan de cette variable. L'intégrale V de (6) est donc une fonc- 
tion entière de k, et il en est alors de même de l'intégrale U de (6), comme nous 
voulions l'établir. Les mêmes raisonnements sont valables pour les séries 
formées avec les dérivées partielles du premier ordre. 

3. Passons au cas de l'équation non linéaire 

8^-^C*'^'^'a^' a^;- 

Nous abrégerons l'exposition, sans diminuer la généralité, en supposant que 
z=0 pour x = 0, et aussi pour «/ = 0. Il suffit évidemment pour cela de rem- 
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placer z par z -f [4» (œ) + 4' (2/) — 4> (0)] • Ceci posé, nous admettons que la 

fonction 

r(a;, y, z, u, v) 

est continue quand {x, y) est dans R, quand z varie entre — a et + a, et que m 
et V varient entre — S et + 5. De plus, pour œ, ?/, z, u et v dans ces intervalles, 
on a 

|F(x, ?/, s', < -w') — F(a;, 2/,s, M, v)|<^i|2;' — s| + ^gjw' — t«| + A;3|w' — ■», 

les h étant des constantes positives Soit enfin M le maximum de la valeur 
absolue de F dans la région où cette fonction est définie. 
On considère les équations successives 

^'"1 -:F{x,y,o,o,o), 



dxdy 

dx dy ^ V^"^ ' '"^ ' "1 ' 3a» ' ^.y 



* (^«^..V, ^1. 3^ , 9-;, 






les z s'annulant tous pour x = aro, quel que soit y, et pour y^yo quel que soit «. 
On sera assuré que 

dx ' dy 

restent compris dans les limites indiquées, si (ce, y) est à l'intérieur d'un rect- 
angle compris dans R, ayant pour sommet (xo, yo), et dont les côtés p et p' satis- 
font aux inégalités 

Mpp'<a, Mp<5, Mp'<&. 

Nous supposerons d'ailleurs que p et p' sont au plus égaux aux côtés a et /3 du 
rectangle R. 

Dans ces conditions la série 

% ~r ^2 "T • • • • ~r 2,j + « • • . 

représentera l'intégrale cherchée. On est, en effet, ramené immédiatement au 
cas de l'équation linéaire, en considérant les équations 

^ = F(.,y,0,0,0). 
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et en leur substituant les équations linéaires 



dxdy ^ "-^ ' ' dx ' ' dy ' 

La convergence de la suite déduite de ces dernières équations entraîne 
immédiatement la convergence de la série des z dans le rectangle (p, p'), et le 
problème est par suite résolu. L'intégrale est déterminée dans le rectangle (p, p'). 

4. Il est remarquable que, dans la question précédente, les limites trouvées 
pour p et p' ne dépendent pas des constantes h. Il faut cependant que l'on soit 
assuré de l'existence de ces constantes pour que le raisonnement soit valable. Un 
cas intéressant est celui où la fonction 

F (a;, y, z, u, v) 

serait déterminée et continue pour toute valeur réelle de z, u et v [le point (x, y) 
étant dans le rectangle R] et où cette fonction aurait des dérivées premières 

ap aF 8F 

dz ' du ' dv 

restant en valeur absolue moindre qu'un nombre fixe dans les mêmes conditions. 

Nous n'aurons pas alors à nous préoccuper des inégalités (7), puisque la 
fonction F est déterminée pour toute valeur de z, u, v; par suite, dans ce cas, 
la série représentant l'intégrale convergera dans R . 

Ainsi, par exemple, l'équation 

dh dz . -, dz . 

'=-a^ V o -^- +csms, 



dx dy dx dy 

où a, &, G sont des fonctions continues de a? et ^Z dans le rectangle R, admettra 
ce rectangle même comme champ de convergence pour la série donnée par les 
approximations successives. On sait que l'équation 



d\ 



sin z 



dxdy 

se rencontre dans la théorie des surfaces à courbure constante, et, dans ses lepons 
de Géométrie différentielle, M. Blanchi s'est servi des approximations successives 
appliquées à cette équation pour traiter un intéressant problème de Géométrie. 

5. Bien d'autres problèmes concernant les équations précédentes pourraient 
être traités par une autre voie analogue. Pour indiquer au moins un nouvel 
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exemple reprenons l'équation linéaire 

ox oy ox oy 
et construisons sur Oœ et Oy (Fig. 89) le carré OABC de côtés OA = 00 = a, 

FiG. 89. 



y 
c 








/ 


p 


/ 


/ 

M 


B 







a 


1 A X 



que nous désignerons par R. On se donne la valeur d'une intégrale z sur OA et 

sur OB ; on aura ainsi 

a = /(x) pour ,7 = 0, 

z = 4)(œ) pour î/ = a;; 

/(x) et 4)(x) sont deux fonctions continues ainsi que leurs dérivées du premier 
ordre j elles sont définies dea; = 0àa3 = a, et l'on a, bien entendu, /(O) = 4) (0) . 

L'intégrale de l'équation 

3^2 



dxdy 



0, 



satisfaisant à ces conditions initiales, sera évidemment 

z=/{x) + f(y)-/{y). 

Ensuite, en désignant par P (a;, y) une fonction donnée de x et de y dans le carré 
R, l'intégrale de l'équation 

s'annulant sur OA et sur OB , sera 

le champ d'intégration est le rectangle MPQR, en désignant par M le point 

(ie, y)- 

Formons alors, comme précédemment, le système 

dx dy 

dhi 9zi 



dx dy dx dy ^ ' 

dK --, gg»-i _. ^ 9a„_i , 
X oy dx dy 
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On intègre la première avec les conditions 

2!i=/(a') pour 2/=0 et Si = 4)(a3) pour y=-x, 

et pour n supérieur kun, on prend 

2» = pour 2/ = et z„ = pour y=-x. 

Des considérations analogues à celles que nous avons employées ci-dessus permet- 
tent aisément d'établir que la série 

^1 ~r ^2 ~r • • • • r ^n ~f" • • • • 

converge uniformément dans R, et représente l'intégrale cherchée. 

6. Comme exemple d'équations d'ordre supérieur au second, pour lesquelles 
s'appliquent sans difficultés les méthodes précédentes, je citerai les équations 
suivantes étudiées à ce point de vue par M. Delassus dans un des chapitres de 
son intéressante thèse {voir aussi Comptes rendus, 1893). Ce sont les équations 
d'ordre n de la forme 

avec les conditions suivantes : 

* = , 1 , . . . . « , , _i_ „\ 

fcO.l s (i. + « = «. i^+0) 

et en supposant A^^ = 1 . 

7. Je voudrais maintenant considérer des équations pour lesquelles on ne 
puisse appliquer la méthode précédente d'approximations. Il n'est pas difficile 
de trouver de tels exemples, nous n'avons qu'à prendre l'équation du premier 
ordre 

dz , . dz ,„, 

^-a{x,y)^. (8) 

Supposons qu'on veuille trouver l'intégrale de cette équation se réduisant, 
pour a; = «o , à une fonction donnée F {y). On peut former les équations suivantes 

^1 = 
dx ' 

dx ''^''^y^dy^ 

dx ""^"^^yi dy ' 
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% prenant pour x = Xq la valeur F {y), et les autres z s'annulant identiquement 
pour cette valeur de x. Mais on voit que l'on ne pourra former les fonctions 
23, ... , z„, . . . que si F (jy) et a (aj, y) ont des dérivées partielles de tout ordre 
par rapport à «/ > et la convergence du développement ne peut être établie que si 
l'on suppose que F {y) et a (a; , y) sont des fonctions analytiques. Il semble donc 
qu'on ne puisse établir V existence des intégrales de V équation (8) qu'en admettant que 
a{x, y) est analytique. Quoique la question n'ait qu'un intérêt théorique, elle 
vaut peut-être la peine d'être examinée. 

Reprenons d'abord, à cet effet, l'étude de l'équation différentielle ordinaire 

i-=-^('^'^)' 

et plapons-nous dans les hypothèses du n° 1. Nous avons trouvé une intégrale 
prenant pour œ = Xo la valeur yo > soit 

y = ¥{x, xo, yo), 

en mettant en évidence toutes les quantités dont dépend F. La fonction F est une 
fonction continue de ce , Xq et y^ ; elle a une dérivée première par rapport à x , 
mais toute la difficulté de la question qui nous occupe est de savoir si cette fonc- 
tion a une dérivée partielle du premier ordre par rapport à 2/0 • Or les approxi- 
mations successives nous conduisent à la suite de fonctions 

se calculant de proche en proche au moyen des formules 

2/1 = 2/0+ J^/i^ ♦ 2/0) <^x , 



2/»= 2/0 -^J^/i^' Vn-i) dx, 

et F (œ , cco , ^o) ^st la limite de 2/» • On peut calculer de proche en proche les 

dérivées partielles 

^ ?^ ... ^ 
92/0 ' 32/o ' 92/o ' 

si l'on admet seulement que /(as, y) a une dérivée partielle du premier ordre par 
rapport ^y . Il est donc bien vraisemblable que F aura une dérivée partielle du 
13 
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premier ordre par rapport à 2/0 • Nous le démontrerons élégamment sans calculs 
en rattachant la question à un problème traité plus haut ; on va supposer que 
f{x, y) a des dérivées partielles des deux premiers ordres par rapport à y. 

fil suffirait même d'admettre que la dérivée -=^ , sans avoir de dérivée par rap- 
port à y, jouit de la propriété dont jouissait la fonction appelée /(ce, y) au n° 1.1 

Je considérerai, dans ce qui va suivre, a^o comme une constante numérique 
et ^ désignera une seconde quantité numérique. J'envisage l'équation aux 
dérivées partielles 

9^2/ _ 9/ ^y 



dx dyo dy dyo ' 



(9) 



définissant une fonction y des deux variables x et y^. D'après ce qui a été vu au 
n° 3, nous pouvons l'intégrer en prenant les conditions initiales suivantes : 

2/ = «/o pour !r = a-o, 

2/ = F(a;, »«, /3) pour 2/0 = /?- 

L'intégrale y de l'équation (9) sera alors complètement déterminée. Or, on 
déduit de cette équation 

P (a;) ne dépendant pas de y^. Or, pour y„=:l3 , on a 

puisque, pour yQ = (3, y est l'intégrale de l'équation -^=/(a;, y), qui prend, 

pour x = Xo, la valeur (3 . La fonction P (x) est donc identiquement nulle, et 
l'intégrale 

2/(a;, 2/0), 

que nous venons d'obtenir, est l'intégrale de l'équation -^=f(x, y) prenant 

pour x = Xo la valeur y^; elle est identique à la fonction F (œ, x^, y^, et celle-ci 
a, par suite, une dérivée du premier ordre par rapport à y^ . 

On démontrera d'une manière analogue que F (», x^, y^ a une dérivée du 
premier ordre par rapport à ajQ. On regardera y^ comme une constante numérique, 
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(10) 



et soit a une seconde quantité numérique. Nous formons l'équation 

9V — 9/ 9y 
dx dxo dy dx^ ' 

en l'intégrant avec les conditions initiales 

y = yo pour x = Xo, 

y = F{x, a,yo) pour Xo = a, 

ce qui correspond au cas étudié (n° 5). On déduit de l'équation (10) 

^=f{x,y) + Qix). 

Q (x) ne dépendant pas de Xq ; on voit que Q (ce) est nul, en faisant dans cette 
relation £Co = a , et l'on termine comme plus haut. 

Ainsi, la fonction F (x, Xo, 2/0) ^ des dérivées partielles du premier ordre par 
rapport h x^ et h yg. Or, la relation 

y = ¥{x, Xq, yo) 

peut manifestement s'écrire 

2/0 = F (xo, X, y), 

puisque l'intégrale qui, pour la valeur x de la variable, prend la valeur y aura 
en Xq la valeur y^. D'après ce qui précède, F(a;o, x, y) est une fonction continue 
de ce et y, et elle a des dérivées partielles du premier ordre elles-mêmes con- 
tinues. Désignons cette fonction par 

F (ce, y), 
en n'écrivant plus la constante Xq : nous aurons l'intégrale générale de l'équation 

^ =/(»'. y) sous la forme 

F (x , 2/) = const., 

et F satisfera à l'équation aux dérivées partielles 

aF , .. x9F_. 

Nbv^ avons donc établi l'existence d'tme intégrale de cette équation, et par suite de 
toutes les intégrales, en supposant seulement que f{x, y) est continue et a des dérivées 
partielles des deux premiers ordres par rapport h y . 
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Ainsi, comme application, on peut établir l'existence d'un facteur intégrant 
pour l'expression 

dy + F{x, y) dx, 

en supposant seulement que la fonction P (a; , y) est continue et a des dérivées 
partielles des deux premiers ordres par rapport à y. 

8. Tout ce que nous venons de dire subsistera évidemment si les fonctions 
considérées, au lieu d'être réelles, sont des fonctions complexes des deux variables 
réelles x et y. En particulier, le théorème relatif au facteur intégrant qui vient 
d'être énoncé s'applique aussi bien si l'on a 

P(a;, y)=p{x,y) + iq{x, y), 

petq étant des fonctions réelles de x et y, jouissant des propriétés indiquées. 

Une application, qui offre quelque intérêt, se présente immédiatement. C'est 

une proposition élémentaire, qu'une surface analytique peut être représentée sur 

un plan, de manière qu'il y ait conservation des angles : on a ainsi une carte 

géographique de la surface. La démonstration bien connue de ce théorème 

s'appuie essentiellement sur ce que la surface est analytique ; elle revient à la 

recherche d'un facteur intégrant. Avec l'extension donnée à cette dernière 

recherche, nous n'avons plus besoin d'admettre que la surface est analytique. 

Soit une surface pour laquelle le carré de l'élément linéaire se mette sous la 

forme 

ds^ = Edx^+ 2¥dx dy + G:dy\ 

On peut, d'après ce que nous venons de dire, démontrer la possibilité de faire la 
carte de cette surface sur un plan, si les trois coefficients E, F, G sont des fonc- 
tions continues de x et y, ayant des dérivées partielles des deux premiers ordres 
par rapport à ?/. Il suffit même de supposer que les trois dérivées du premier 
ordre 

aE aF 30 

dy ' dy ' dy ' 

jouissent de la propriété admise pour la fonction / (ce , y) au n° 1. 

Ces conditions, un peu dissymétriques, sont suffisantes ; elles ne sont sans 
doute pas toutes nécessaires, mais, pour s'en affranchir, il faudrait trouver un 
autre mode de démonstration pour l'existence du facteur intégrant. 



